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Abstract 

This  is  a  primer  on  extended  Gaussian  Images.  Extended  Gaussian  Images  are 
useful  for  representing  the  shapes  of  surfaces.  They  can  be  computed  easily  from: 

1.  Needle  maps  obtained  using  photometric  stereo,  or 

2.  Depth  maps  generated  by  ranging  devices  or  stereo. 

Importantly,  they  can  also  be  determined  simply  from  geometric  models  of  the 
objects.  Extended  Gaussian  images  can  be  of  use  in  at  least  two  of  the  tasks  facing 
a  machine  vision  system: 

1.  Recognition,  and 

2.  Determining  the  attitude  in  space  of  an  object. 

Here,  the  extended  Gaussian  image  is  defined  and  some  of  its  properties  discussed. 
An  elaboration  for  non-convex  objects  is  presented  and  several  examples  are  shown. 
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1.  Introduction 


In  order  to  recognize  an  object  and  to  determine  its  attitude  in  space,  it  is 
necessary  to  have  a  way  of  representing  the  shape  of  its  surface.  Giving  the  distance 
to  the  surface  along  parallel  rays  on  a  regularly  spaced  grid  provides  one  way  of 
doing  this.  This  simple  representation  is  called  a  depth  map.  A  range  finder  produces 
surface  descriptions  in  this  form.  Unfortunately,  depth  maps  do  not  transform  in  a 
simple  way  when  the  object  rotates  (For  one  thing,  interpolation  must  be  used  to 
get  a  new  depth  map  on  a  regularly  spaced  grid).  Alternatively,  surface  orientation 
might  be  given  for  points  on  the  surface  on  some  regular  sampling  grid.  This 
grid  may  conveniently  correspond  to  the  picture  cells  in  an  image.  Such  a  simple 
representation  is  called  a  needle  map  (Figure  1)  [13].  Photometric  stereo  produces 
surface  descriptions  in  this  form  [7,  12,  15,  28,  32,  33 j .  This  representation  too  is 
not  directly  helpful  when  it  comes  to  comparing  surfaces  of  objects  that  may  be 
rotated  relative  to  one  another  (Both  depth  maps  and  needle  maps  depend  on  the 
position  of  the  object  as  well  as  its  attitude). 

The  extended  Gaussian  image,  on  the  other  hand,  does  make  it  easy  to  deal 
with  the  varying  attitude  of  an  object  in  space  [2,  3,  9,  13,  16,  29].  For  one  thing, 
it  is  insensitive  to  the  position  of  the  object.  Some  information  appears  to  be 
discarded  in  the  formation  of  the  extended  Gaussian  image,  however.  Curiously,  in 
the  case  of  convex  objects,  the  representation  is  nevertheless  unique.  That  is,  no 
two  convex  objects  have  the  same  extended  Gaussian  image. 

This  representation  of  the  shape  of  the  surface  of  an  object  has  proven  most 
useful  in  work  on  automatic  bin  picking  [4,  17,  18].  We  start  our  discussion 
with  objects  having  planar  faces.  Later  we  consider  smoothly  curved  objects. 
Methods  for  computing  discrete  approximations  of  extended  Gaussian  images, 
called  orientation  histograms,  arc  presented  too.  Orientation  histograms  can  be 
computed  from  experimental  data  or  -mathematical  descriptions  of  the  objects. 
Sections  marked  with  an  asterisk  may  be  omitted  on  first  reading  or  if  your  interest 
in  the  mathematical  details  is  limited. 

2.  Discrete  Case:  Convex  Polyhedra 

Minkowski  showed  in  1897  that  a  convex  polyhedron  is  fully  specified  (up  to 
translation)  by  the  area  and  orientation  of  its  faces  [21,  23,  26].  We  can  represent 
area  and  orientation  of  the  faces  conveniently  by  point  masses  on  a  sphere.  Imagine 
moving  the  unit  surface  normal  of  each  face  so  that  its  tail  is  at  the  center  of  a  unit 
sphere.  The  head  of  the  unit  normal  then  lies  on  the  surface  of  the  unit  sphere.  This 
sphere  is  called  the  Gaussian  sphere  and  each  point  on  it  corresponds  to  a  particular 
surface  orientation.  The  extended  Gaussian  image  of  the  polyhedron  is  obtained 
by  placing  a  mass  at  each  point  equal  to  the  surface  area  of  the  corresponding  face 
(Figure  2). 

It  seems  at  first,  as  if  some  information  is  lost  in  this  mapping,  since  the 
position  of  the  surface  normals  is  discarded.  Viewed  another  way,  no  note  is 
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Figure  1.  A  needle  map  shows  unit  surface  normals  at  points  on  the  surface  on 
a  regular  grid.  Normals  which  point  towards  the  viewer  will  be  seen  as  dots,  while 
tilted  surface  patches  give  rise  to  normals  which  are  shown  as  lines  pointing  in  the 
direction  of  steepest  descent. 


made  of  the  shape  of  the  faces  or  their  adjacency  relationships.  It  can  nevertheless 
be  shown  that  (up  to  translation)  the  extended  Gaussian  image  uniquely  defines 
a  convex  polyhedron  [29].  An  iterative  algorithm  has  recently  been  invented  for 
recovering  a  convex  polyhedron  from  its  extended  Gaussian  image  [22]. 

2.1.  Properties  of  the  Extended  Gaussian  Image 

The  extended  Gaussian  image  is  not  affected  by  translation  of  the  object. 
Rotation  of  the  object  induces  an  equal  rotation  of  the  extended  Gaussian  image, 
since  the  unit  surface  normals  rotate  with  the  object. 

Mass  distributions  which  lie  entirely  within  one  hemisphere,  that  is,  are  zero 
in  the  complementary  hemisphere,  do  not  correspond  to  closed  objects.  As  we  shall 
see,  the  center  of  mass  of  an  extended  Gaussian  image  has  to  lie  at  the  origin.  This 
is  clearly  not  possible  if  a  whole  hemisphere  is  empty.  Also,  a  mass  distribution 
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Figure  2.  The  extended  Gaussian  image  of  a  polyhedron  can  be  thought  off 
as  a  collection  of  point  masses  on  the  Gaussian  sphere.  Each  mass  is  proportional 
to  the  area  of  the  corresponding  face.  Point  masses  on  the  visible  hemisphere  are 
shown  as  solid  marks,  the  others  as  open  marks.  The  center  of  mass  (shown  as  the 
symbol  0)  must  be  at  the  center  of  the  sphere  if  the  polyhedron  is  a  closed  object. 


which  is  non-zero  only  on  a  great  circle  of  the  sphere  corresponds  to  the  the  limit 
of  a  sequence  of  cylindrical  objects  of  increasing  length  and  decreasing  diameter 
(Figure  3).  We  will  exclude  such  pathological  cases  and  confine  our  attention  to 
closed,  bounded  objects  [23,  29]. 

Some  properties  of  the  extended  Gaussian  image  are  important:  Firstly,  the 
total  mass  of  the  extended  Gaussian  image  is  obviously  just  equal  to  the  total 
surface  area  of  the  polyhedron.  If  the  polyhedron  is  closed,  it  will  have  the  same 
projected  area  when  viewed  from  any  pair  of  opposite  directions.  This  allows  us  to 
compute  the  location  of  the  center  of  mass  of  the  extended  Gaussian  image. 

Imagine  viewing  a  convex  polyhedron  from  a  great  distance.  Let  the  direction 
from  the  object  towards  the  viewer  be  given  by  the  unit  vector  v.  A  face,  with  unit 
norma!  s,,  will  be  visible  only  if  s,  •  v  >  0.  Suppose  that  the  surface  area  of  this 
face  is  O,.  Due  to  foreshortening  it  will  appear  only  as  large  as  would  a  face  of  area 

(S,-v)Ot, 

normal  to  v  (Figure  4).  The  total  apparent  area  of  the  visible  surface  is 

A(v)=  £  (Si  •  v)Oi, 

{»  I  S,  V>0} 
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when  viewed  from  the  direction  v.  The  total  apparent  area  of  the  visible  surface 
when  viewed  from  the  opposite  direction  is 

A(— v)  =  •  E  (»i  •  v)  Oi, 

{i  |s,-v<0> 

This  should  be  the  same,  that  is,  A(v)  =  A( — v).  Consequently, 

E  (■»  •  *)  Oi  =  E  s,  o,  •  v  =  o, 

all  t  Lail  » 

where  the  sum  now  is  over  ail  faces  of  the  object.  This  holds  true  for  all  view 
vectors,  v,  so  we  must  have 

E  5«  Oi  =  o. 

all  i 

That  is,  the  center  of  mass  of  the  extended  Gaussian  image  is  at  the  origin. 

An  equivalent  representation,  called  a  spike  model,  is  a  collection  of  vectors 
each  of  which  is  parallel  to  one  of  the  surface  normals  and  of  length  equal  to  the 
area  of  the  corresponding  face.  The  result  regarding  the  center  of  mass  is  equivalent 
to  the  statement  that  these  vectors  must  form  a  closed  chain  when  placed  end  to 
end  (Figure  5). 


m 


Figure  -1.  A  surface  element  appears  smaller  because  of  foreshortening.  The 
apparant  area  is  the  true  area  times  the  cosine  of  the  angle  between  the  surface 
normal  and  the  vector  pointing  towards  the  viewer. 


Figure  5.  Vectors  parallel  to  the  normals  of  the  faces  of  a  polyhedron,  and 
of  length  equal  to  the  areas  of  the  corresponding  faces,  form  a  closed  chain  when 
placed  end  to  end. 


2.2.  Reconstruction  of  a  Tetrahedron  (*) 

Faces  that  share  a  common  edge  are  said  to  be  adjacent.  The  masses  on  the 
Gaussian  sphere  corresponding  to  two  adjacent  faces  are  not  always  each  others 
closest  neighbors.  Recovering  a  polyhedron  from  its  extended  Gaussian  image  is 
not  easy  in  the  general  case  [22],  because  it  is  hard  to  determine  which  faces  are 
adjacent.  Pending  the  actual  offsets  of  each  of  the  faces  from  the  center  of  mass  of 
the  polyhedron  is  not  as  hard. 


The  structure  of  a  tetrahedron,  however,  is  very  simple:  Every  face  is  adjacent 
to  the  other  three.  The  shape  of  the  tetrahedron  is  completely  determined  by  the 
surface  normals  of  the  four  faces,  only  the  size  of  the  tetrahedron  remaining  to  be 
determined.  In  other  words:  There  is  only  one  degree  of  freedom  left.  Another  way 
to  look  at  it  is  to  note  that  the  four  faces  must  have  areas  that  place  the  center 
of  mass  of  the  extended  Gaussian  image  at  the  origin,  as  we  have  just  seen.  This 
condition  places  three  constraints  on  the  four  parameters. 

Let  the  given  unit  surface  normals  be  a,  b,  c,  and  d,  and  the  areas  of  the 
corresponding  faces,  A,  B,  C,  and  D  (Figure  6).  We  have  to  determine  the  distances, 
a,  6,  c,  and  d,  of  these  faces  from  the  center  of  mass  of  the  tetrahedron.  From  these 
distances  we  can,  if  desired,  compute  the  positions  of  the  vertices  A,  B,  C,  and  D, 
simply  by  intersecting  three  of  the  planes  at  a  time.  The  notation  here  is  that  the 
face  opposite  vertex  A  has  area  A  and  unit  surface  normal  a,  and  so  on. 

The  perpendicular  distance  of  the  center  of  area  of  a  triangle  from  one  of 
the  sides  is  one  third  the  perpendicular  distance  of  the  vertex  opposite  that  side. 
Similarly,  in  a  tetrahedron,  the  distance  from  the  center  of  mass  to  a  particular  face 
is  equal  to  one  quarter  of  the  distance  of  the  vertex  opposite  that  face.  We  start  by 
Finding  a  formula  for  the  distance  of  the  face  with  area  D,  say,  from  the  opposite 
vertex  d.  The  desired  distance,  d,  will  be  just  a  quarter  of  the  result  obtained  in 
this  fashion.  The  remaining  three  distances,  a,  b,  and  c,  can  then  be  computed 
using  formulae  obtained  by  cyclical  permutation  of  the  variables. 
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I  he  position  of  the  reconstructed  tetrahedron  is  abritrary,  since  the  extended 
Gaussian  image  is  insensitive  to  translation.  To  make  the  result  unique,  we  might 
place  the  center  of  mass  at  the  origin.  To  reduce  the  size  of  the  expressions  to  be 
manipulated  here,  however,  it  is  convenient  to  move  the  tetrahedron  so  that  one 
vertex,  D,  say,  is  at  the  origin.  The  distances  of  the  faces  from  the  center  of  mass 
are  obviously  not  affected  by  this. 

Suppose  for  now  that  we  know  the  locations  of  the  vertices  A,  B,  and  C  relative 
to  I).  We  can  then  compute  the  directions  of  the  six  edges  of  the  tetrahedron, 
by  taking  all  of  the  distinct  pairwise  differences  of  the  four  vertex  positions.  Four 
surface  normals  can  then  be  found  by  taking  cross-products  of  these  edge-direction 
vectors.  We  actually  need  only  four  of  the  edge  vectors  forming  a  closed  circuit  to 
do  this.  The  results  can  then  be  normalized  to  obtain  unit  surface  normals, 

a  B  XC  .  _CXA  . _  A  X  B 
B  XC|!  “  |C  XA]’  ~~  j  A~X  B]  ’ 

and, 


*  (A  -  C)  X  (B  -  A)  AXB  +  BXC  +  CXA 
'  1(A  —  C)  X  (B  -  A)|  “  |A  X  B  +  B  X  C  +  C  X  A| ' 

Now  the  perpendicular  distance  of  the  plane  with  area  D  from  the  origin  can  be 
found  by  taking  the  dot-product  of  any  of  the  three  vertices,  A,  B,  and  C  with  the 
unit  normal  d.  Thus, 


4d  =  dA  =  dB  =  dC  = 


[ABC] 

|  A  X  B  +  BxC  +  CxAf 


The  area  of  the  facet  opposite  the  origin  is  also  easy  to  compute, 


D  —  |  [(A  —  C)  X  (B-A)|  =  ||A  X  B  +  B  X  C-f  C  X  A|. 

Our  task  is  to  express  the  oITset  d  in  terms  of  the  area  D  and  the  given  unit  surface 
normals.  The  two  formulae  above  do  not  allow  us  to  do  that  directly,  because  we 
do  not  know  what  the  value  of  [ABC]  is.  This  quantity,  by  the  way,  is  six  times 
the  volume,  V ,  of  the  tetrahedron,  or, 


v-^M0=5|abc), 

We  proceed  by  considering  the  four  distinct  triple  products  of  the  four  unit  surface 
normals.  First  of  all, 
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[ABC]2 


[A  X  B||B  x  C|jC  X  A|  ’ 
since  [(x  X  y)(y  X  z)  (z  X  x)]  =  [xyz]2.  Then,  by  similar  reasoning, 

[ABC]2 


[abd]  = 


|B  X  C||C  X  A|  |  A  XB  +  BXC  +  CXAI’ 
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since  [x  y  (x  +  y  -f  i)]  =  [xyz].  Formulae  for  [bed]  and  [cad]  can  be  found  by 
cyclical  permutation  of  the  variables. 


Multiplying  the  three  formulae  found  this  way  together  we  get, 


[abd][bcd][cad] 
and  so 


[ABC]6 

(A  X  B)2(B  X  C)2(C  X  A)2 1 A  X  B-f  B  x  C  +  C  X  A|3’ 


[abd][bcd][cad] 
[a  be]2 

So  that  finally, 


[ABC]2 

|AXB  +  BXC  +  C  X  A| 


(4<0’(2fl). 


y 1  (2D)[a  b  d] [b  c  d] [c  a  d] 

—  [a  be] 

The  other  distances,  a,  b,  and  c,  can  be  computed  using  similar  fo  ae  obtained 
by  cyclical  permutation  of  the  variables. 


3.  Continuous  Case:  Smoothly  Curved  Objects 

The  ideas  presented  in  the  previous  chapter  can  be  extended  to  apply  to 
smoothly  curved  surfaces. 

3.1.  Gaussian  Image 

One  can  associate  a  point  on  the  Gaussian  sphere  with  a  given  point  on  a 
surface  by  finding  the  point  on  the  sphere  which  has  the  same  surface  normal 
(Figure  7)  [6,  11,  23].  Thus  it  is  possible  to  map  information  associated  with  points 
on  the  surface  onto  points  on  the  Gaussian  sphere.  There  are  no  two  points  with 
the  same  surface  normal  in  the  case  of  a  convex  object  which  has  positive  Gaussian 
curvature  everywhere.  The  mapping  from  the  object  to  the  Gaussian  sphere  in  this 
case  is  invertible:  Corresponding  to  each  point  on  the  Gaussian  sphere  there  is  a 
unique  point  on  the  surface  (If  the  convex  surface  has  patches  with  zero  Gaussian 
curvature,  curves  or  even  areas  on  it  may  correspond  to  a  single  point  on  the 
Gaussian  sphere). 

One  useful  property  of  the  Gaussian  image  is  that  it  rotates  with  the  object. 
Consider  two  parallel  surface  normals,  one  on  the  object  and  the  other  on  the 
Gaussian  sphere.  The  two  normals  will  remain  parallel  if  the  object  and  the  Gaussian 
sphere  are  rotated  in  the  same  fashion.  A  rotation  of  the  object  thus  corresponds 
to  an  equal  rotation  of  the  Gaussian  sphere. 

3.2.  Gaussian  Curvature 

Consider  a  small  patch  SO  on  the  object.  Each  point  in  this  patch  corresponds 
to  a  particular  point  on  the  Gaussian  sphere.  The  patch  SO  on  the  object  maps 
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Figure  7.  The  Gaussian  image  of  an  object  is  obtained  by  associating  with 
each  point  on  us  surface  the  point  on  the  Gaussian  sphere  which  has  the  same 
surface  orientation.  The  mapping  is  invertible  if  the  object  has  positive  Gaussian 
c  u  r va t  u  re  e ve  r y w h  c re . 


into  a  patch,  6S  say.  on  the  Gaussian  sphere  (Figure  8).  If  the  surface  is  strongly 
curved,  the  normals  of  points  in  the  patch  will  point  into  a  wide  fan  of  directions. 
The  corresponding  points  on  the  Gaussian  sphere  will  be  spread  out.  Conversely,  if 
the  surface  is  planar,  the  surface  normals  are  parallel  and  map  into  a  single  point. 

Those  considerations  suggest  a  suitable  definition  of  curvature.  The  Gaussian 
curvature  is  defined  to  be  equal  to  the  limit  of  the  ratio  of  the  two  areas  as  they 
tend  to  zero.  That  is, 


K 


,.  6S  dS 
hm  ■—  =  — . 
<50-0  <fO  dO 


From  this  differential  relationship  we  can  obtain  two  useful  integrals.  Consider  first 
integrating  K  over  a  finite  patch  O  on  the  object: 


L!  Kdo= U ds=s- 

where  S  is  the  area  of  the  corresponding  patch  on  the  Gaussian  sphere.  The 
expression  on  the  left  is  called  the  integral  curvature.  This  relationship  allows  one 
to  deal  with  surfaces  which  have  discontinuities  in  surface  normal. 


Now  consider  instead  integrating  \/K  over  a  patch  S  on  the  Gaussian  sphere: 
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Figure  8.  A  patch  on  the  object  maps  into  a  patch  on  the  Gaussian  sphere. 
The  Gaussian  curvature  is  the  limit  of  the  ratio  of  the  area  of  the  patch  on  the 
Gaussian  sphere  to  the  area  of  the  patch  on  the  object  as  these  become  smaller  and 
smaller. 


fsJl/KdS  =  JoJd°  =  °, 

where  O  is  the  area  of  the  corresponding  patch  on  the  object.  This  relationship 
suggests  the  use  of  the  inverse  of  the  Gaussian  curvature  in  the  definition  of  the 
extended  Gaussian  image  of  a  smoothly  curved  object,  as  we  shall  see.  It  also 
shows,  by  the  way,  that  the  integral  of  1  /K  over  the  Gaussian  sphere  equals  the 
total  area  of  the  object. 

3.3.  Alternate  Definition  of  Gaussian  Curvature 

Consider  a  plane  which  includes  the  surface  normal  at  some  point  on  a  smooth 
surface.  The  surface  cuts  this  plane  along  a  curve  called  a  normal  section  (Figure 
9)  [6,  11,  20,  24,  26].  Let,  the  curvature  of  the  normal  section  be  denoted  by 
k.\-.  Consider  the  one-parameter  family  of  planes  containing  the  surface  normal. 
Suppose  that  0  is  the  angle  between  a  particular  plane  and  a  given  reference  plane. 
Then  varies  with  0  in  a  periodic  fashion.  In  fact,  if  we  measure  0  from  the  plane 
that  gives  maximum  curvature,  then  it  can  be  shown  that 

x^(9)  =  cos2  9  -f-  K2  sin2  9, 

where  k\  is  the  maximum  and  ^2  is  the  minimum  curvature.  These  two  values 
of  k  are  called  the  principal  curvatures.  The  corresponding  planes  are  called 
the  principal  planes.  The  two  principal  planes  are  orthogonal,  provided  that  the 
principal  curvatures  are  distinct  (Figure  9). 
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Figure  9.  Norma!  sections  of  the  surface  are  made  with  planes  which  include 
the  surface  normal.  The  planes  corresponding  to  the  largest  and  smallest  value  of 
curvature  are  referred  to  as  the  principal  planes.  The  Gaussian  curvature  is  equal 
to  the  product  of  the  largest  and  smallest  values  of  curvature. 


It  turns  out  that 

K  =  Ki«2 

is  equal  to  the  Gaussian  curvature  introduced  earlier.  This  is  clearly  zero  for  a 
plane.  It  is  equal  to  1  /  R2  for  a  spherical  surface  of  radius  R,  since  the  curvature 
of  any  normal  section  is  1  /R. 

A  ruled  surface  is  one  which  can  be  generated  by  sweeping  a  line  through 
space.  A  hyperboloid  provides  one  example  of  such  a  surface.  Developable  surfaces 
are  special  cases  of  ruled  surfaces  [6,  11,  24].  Cylindrical  and  conical  surfaces  are 
examples  of  developable  surfaces  (Figure  10).  On  a  developable  surface  at  least  one 
of  the  two  principal  curvatures  is  zero  at  all  points.  Consequently  the  Gaussian 
curvature  is  zero  everywhere  too. 

3.4.  The  Extended  Gaussian  Image 

We  can  define  a  mapping  which  associates  the  inverse  of  the  Gaussian  curvature 
at  a  point  on  the  surface  with  the  corresponding  point  on  the  Gaussian  sphere.  Let 
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Figure  10.  A  conical  surface  is  an  example  of  a  developable  surface.  On  it 
the  Gaussian  curvature  is  everywhere  zero,  because  (at  least)  one  of  the  principal 
curvatures  is  zero. 


k,  and  v  be  parameters  used  to  identify  points  on  the  original  surface.  Similarly,  let 
£  and  r)  be  parameters  used  to  identify  points  on  the  Gaussian  sphere  (These  could 
be  longitude  and  latitude,  for  example).  Then  we  define  the  extended  Gaussian 
image  as 

'  *?)  ==  777  7 » 

K{u,v] 

where  (£,77)  is  the  point  on  the  Gaussian  sphere  with  the  same  normal  as  the 
point  [u,v)  on  the  original  surface.  It  can  be  shown  that  this  mapping  is  unique 
(up  to  translation)  for  convex  objects.  That  is,  there  is  only  one  convex  object 
corresponding  to  a  particular  extended  Gaussian  image  [l,  26,  29].  The  proof  is 
unfortunately  non-constructive  and  no  direct  method  for  recovering  the  object  is 
known. 

3.5.  Properties  of  the  Extended  Gaussian  Image 

The  center  of  mass  of  the  extended  Gaussian  image  of  a  smoothly  curved 
object  is  at  the  origin.  We  show  this  in  a  way  similar  to  that  used  earlier  for 
extended  Gaussian  images  of  polyhedral  objects.  Consider  viewing  a  convex  object 
from  a  great  distance.  Let  the  direction  from  the  object  towards  the  viewer  be 
given  by  the  unit  vector  v.  A  surface  patch  with  unit  normal  s  will  be  visible  only 
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if  s  ■  v  >  0.  Suppose  its  surface  area  is  60  (Figure  4).  Due  to  foreshortening  it  will 
appear  only  as  large  as  would  a  patch  of  area 

(s  •  y)SO, 

normal  to  v.  Let  //(v)  be  the  unit  hemisphere  for  which  s  v  >0.  Then  the  apparent 
area  of  the  visible  surface  is 

Aiy)  =  JHJ  G{h){i-v)dS, 

when  viewed  from  the  direction  v.  The  apparent  area  of  the  visible  surface  when 
viewed  from  the  opposite  direction  is 

*(->)  =  Sn^S  cm- ->)iS. 

This  should  be  the  same,  that  is,  A(v)  =  A( — v).  Consequently, 

fs  J  <?(*)  (s  •  v)  dS  —  ( fs  f  G{ s)  s  ds'j  ■  v  ==  0, 

where  the  integral  now  is  over  the  whole  sphere  5.  This  holds  true  for  all  view 
vectors,  v,  so  we  must  have 

Js  J  G(s)sdS  =  0. 

That  is,  the  center  of  mass  of  the  extended  Gaussian  image  is  at  the  origin  (This, 
by  the  way,  is  not  a  very  helpful  constraint  in  practice,  since  one  usually  only  sees 
one  side  of  the  object). 

Another  property  of  the  extended  Gaussian  image  is  also  easily  demonstrated. 
The  total  mass  of  the  extended  Gaussian  image  equals  the  total  surface  area  of  the 
object.  If  one  wishes  to  deal  with  objects  of  the  same  shape  but  differing  size  one 
may  normalize  the  extended  Gaussian  image  by  dividing  by  the  total  mass. 

We  can  think  of  the  extended  Gaussian  image  in  terms  of  mass  density  on 
the  Gaussian  sphere.  It  is  possible  then  to  deal  in  a  consistent  way  with  places 
on  the  surface  where  the  Gaussian  curvature  is  zero,  using  the  integral  of  l/K 
shown  earlier.  A  planar  region,  for  example,  corresponds  to  a  point  mass.  This  in 
turn  corresponds  to  an  impulse  function  on  the  Gaussian  sphere  with  magnitude 
proportional  to  the  area  of  the  planar  region. 

It  is  tempting,  by  the  way,  to  imagine  that  one  can  find  the  attitude  of  an 
object  by  lining  up  the  principal  axes  of  inertia  of  the  observed  extended  Gaussian 
image  and  the  one  computed  from  the  geometric  model  [29].  This  would  be  rather 
straightforward,  requiring  only  the  calculation  of  the  eigenvectors  of  a  three  by 
three  inertia  matrix.  In  practice,  one  typically  has  information  only  about  the 
visible  hemisphere  and  thus  cannot  compute  the  required  first  and  second  moments 
over  the  whole  sphere. 

3.6.  Objects  that  are  not  Convex 

Three  things  happen  when  the  surface  is  non-convex: 
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1.  The  Gaussian  curvature  for  some  points  will  be  negative. 

2.  More  than  one  point  on  the  object  will  contribute  to  a  given  point  on  the 
Gaussian  sphere. 

3.  Parts  of  the  object  may  be  obscured  by  other  parts. 

We  chose  to  extend  the  definition  of  the  extended  Gaussian  image  in  this  case  to 
be  the  sum  of  the  absolute  values  of  the  inverses  of  the  Gaussian  curvature  at  all 
points  having  the  same  surface  orientation, 

This  definition  is  motivated  by  the  method  used  to  compute  the  extended  Gaussian 
image  in  the  discrete  case,  as  we  will  see  later. 

The  above  extension  makes  sense  if  there  are  a  finite,  or  at  most  a  countable, 
number  of  points  on  the  surface  with  the  same  orientation.  At  times,  however,  all 
points  on  a  curve  or  even  an  area  on  the  surface  have  parallel  surface  normals.  In 
this  case  we  may  use, 

C(n)  =  fQ  j  <5(n  —  s)  dO, 

where  n  is  a  unit  vector  on  the  Gaussian  sphere,  while  s  is  a  unit  vector  on  the 
surface  of  the  object.  The  integration  is  over  the  whole  surface  of  the  object  O  and 
6  is  the  unit  impulse  function  defined  on  a  sphere. 

We  can  be  more  specific,  if  we  let  r(u,  v)  be  a  vector  giving  the  point  on  the 
surface  corresponding  to  the  parameters  u  and  v,  then 

tj)  cos  rj  =  J  j  6{£  —  0(u,v),r)  —  d>(u,  u))|ru  X  rv\dudv, 

where  0{u,  v)  and  4>(u,  v )  are  the  latitude  and  longitude  of  the  point  on  the  Gaussian 
sphere  which  has  the  same  orientation  as  the  surface  does  at  the  point  (u,  o).  A 
planar  region  of  area  A  will  thus  contribute  an  impulse  of  weight  A  to  the  extended 
Gaussian  image,  while  a  cylindrical  region  will  give  rise  to  an  impulse  wall  along  a 
great  circle  at  right  angles  to  the  axis  of  the  cylinder.  The  integral  of  the  impulse 
wall  will  be  equal  to  the  area  of  the  cylindrical  region. 

Usually  we  think  of  the  extended  Gaussian  image  as  a  fixed  entity  associated 
with  an  object.  In  the  case  of  non-convex  objects  we  might  want  to  alter  the  definition 
to  include  only  those  parts  of  the  surface  visible  from  a  particular  direction.  This 
would  make  the  (modified)  Gaussian  image  dependent  on  the  view-point.  We  avoid 
this  potential  complication  here. 

3.7.  Examples  of  Extended  Gaussian  Images  (*) 

The  extended  Gaussian  image  of  a  sphere  of  radius  R  is 

G((,r,)  =  R>, 


'  I 
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as  discussed  already. 

Perhaps  slightly  more  interesting  is  the  case  of  an  ellipsoid  with  semi-axes  a, 
b  and  c  lined  up  with  the  coordinate  axes  (Figure  11).  An  equation  for  its  surface 
can  be  written 

O’+GM)’- 

More  useful  for  our  purposes  here  is  a  parametric  form 

x  =acos0cos^, 
y  =6  sin  9  cos  <f>, 
z  =csin^. 

A  normal  at  the  point 

r  =  (a  cos  9  cos  <£,  6  sin  9  cos  <f>,  c  sin  <t>)T 
on  the  surface  is  given  by 

n  =  (be  cos  9  cos  <t>,  ca  sin  9  cos  4>,  ab  sin  ^)r, 
as  will  be  shown  later.  The  Gaussian  curvature  turns  out  to  be  equal  to 

abc  2  fabcl2 

(be  cos  9  cos  <£)2  -f-  (ca  sin  9  cos  <p)2  -f  (a6  sin  ^)2  In3]’ 

where  n2  =  n  •  n. 


If  we  let  £  be  the  longitude  and  t?  be  the  latitude  on  the  Gaussian  sphere,  then 
a  unit  normal  at  the  point  (£,t?)  on  the  sphere  is  given  by  (Figure  12) 

n  =  (cos  £  cos??,  sin  £  cos  rj,  sin  t?)t. 

Now  n  =  nn.  Identifying  terms  in  the  two  expressions  for  surface  normals  at 
corresponding  points  on  the  ellipsoid  and  the  Gaussian  sphere  we  get, 

be  cos  6  cos  <t>  —  n  cos  £  cos  r?, 
ca  sin  9  cos  <i>  =  n  sin  £  cos  r?, 
at  sin  4>  =  n  sin  r?, 

so  that 

n2[(a  cos  £  cos;?)2  -)-  (6cos  £  sin  t?)2  +  (csin  t?)2]  —  (a6c)J, 
and  finally,  substituting  for  n2  in  the  equation  for  K  we  get 


C(ti)=  K-  = 


abc 


[  (a  cos  £  cos  t?)2  (b  sin  £  cos  r?)  2  +  (c  sin  t?)2J  ' 

The  extended  Gaussian  image,  in  this  case,  varies  smoothly  and  has  the  stationary 
values 
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at  the  points  where  r  is  equal  to  (±1,0,  0)r,  (0,  ±1, 0)r,  and  (0,0,  ±l)r,  respectively. 
These  results  can  be  easily  checked  by  sectioning  the  ellipsoid  using  the  xy ,  yz 
and  xx-planes.  The  Gaussian  curvature,  in  this  case,  equals  the  product  of  the 
curvatures  of  the  resulting  ellipses.  One  then  uses  the  fact  that  the  maximum  and 
minimum  curvatures  of  an  ellipse  with  semi-axes  a  and  6  are  a/b 2  and  6/a2. 

Later  we  will  derive  the  extended  Gaussian  image  of  a  torus,  an  object  that  is 
not  convex. 


4.  Discrete  Approximation:  Needle  Maps 

Consider  the  surface  broken  up  into  small  patches  of  equal  area.  Let  there 
be  p  patches  per  unit  area.  Erect  a  surface  normal  on  each  patch.  Consider  the 
polyhedral  object  formed  by  e  intersection  of  the  tangent  planes  perpendicular 
to  these  surface  normals.  It  approximates  the  original  surface.  The  smaller  the 
patches,  the  better  the  approximation. 

The  extended  Gaussian  image  of  the  original  (smoothly  curved)  convex  object 
is  approximated  by  impulses  corresponding  to  the  small  patches.  The  magnitude  of 
each  impulse  is  about  1/p,  corresponding  to  the  area  of  the  patch  it  rests  on  (Figure 
13).  Strongly  curved  areas  will  distribute  their  impulses  over  a  large  region  on  the 
Gaussian  sphere,  while  areas  which  are  nearly  planar  will  have  them  concentrated 
in  a  small  region.  In  fact,  the  number  of  impulses  per  unit  area  on  the  Gaussian 
sphere  approaches  p  times  the  absolute  value  of  the  Gaussian  curvature  as  we  make 
p  larger  and  larger.  This  can  be  shown  using  the  integral  of  1/ K  given  earlier. 

The  tesselation  of  the  surface  can  be  based  on  an  arbitrary  division  into 
triangular  patches  as  long  as  the  magnitude  of  each  impulse  on  the  Gaussian 
sphere  is  made  proportional  to  the  area  of  the  corresponding  patch  on  the  surface. 
Alternatively,  one  can  divide  the  surface  up  according  to  the  division  of  the  image 
into  picture  cells.  In  this  case  one  has  to  take  into  account  that  the  area  occupied  in 
the  image  by  a  given  patch  is  affected  by  foreshortening.  The  actual  surface  area  is 
proportional  to  l/(s,  •  v),  where  s,  is  the  normal  of  the  patch,  while  v  is  the  vector 
pointing  towards  the  viewer  (Figure  4). 

Measurements  of  surface  orientation  from  images  will  not  be  perfect,  since  they 
are  affected  by  the  noise  in  brightness  measurements.  Similarly,  surface  orientations 
obtained  from  range  data  will  be  somewhat  inaccurate.  Consequently  the  impulses 
on  the  Gaussian  sphere  will  be  displaced  ^  little  from  their  true  positions.  The 
expected  density  on  the  Gaussian  sphere  will  nevertheless  tend  to  be  equal  to 
the  inverse  of  the  Gaussian  curvature.  One  cannot,  however,  expect  the  impulses 
corresponding  to  a  planar  surface  to  be  coincident.  Instead,  they  will  tend  to  form 
a  small  cluster.  To  be  precise,  the  effect  of  noise  is  to  smear  out  the  information  on 
the  sphere.  The  extended  Gaussian  image  is  convolved  with  a  smoothing  function 
of  width  proportional  to  the  magnitude  of  the  noise. 
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4.1.  Using  Object  Models 

Extended  Gaussian  images  also  have  to  be  computed  for  surfaces  of  prototypical 
object  models.  In  this  case  it  is  best  to  find  a  convenient  way  to  parameterize  the 
surface  and  break  it  up  into  many  small  patches.  Suppose  the  surface  is  given  in 
terms  of  two  parameters  u  and  v  as  r(u,v).  Then  at  the  point  ( u,v )  we  see  that  r„ 
and  rv  are  two  tangents  (Figure  14).  The  cross  product  of  these  tangents  is  normal 
to  vhe  surface.  The  unit  normal 

.  X  Tj 

n  =  i — - T 

I ru  X  r„| 

allows  us  to  determine  to  which  point  on  the  Gaussian  sphere  this  patch  corresponds. 
Suppose  that  we  divided  the  range  of  u  into  segments  of  size  6u  and  the  range  of 
v  into  segments  of  sLe  6v.  Then  the  area  of  the  patch, 

6 A  =  jru  X  rv| 6u6v, 

can  be  used  to  determine  what  contribution  it  makes  to  the  corresponding  place  on 
the  Gaussian  sphere.  Note  that  we  do  not  have  to  explicitly  compute  the  Gaussian 
curvature  or  take  second  partial  derivatives. 


Figure  14.  A  surface  normal  can  he  computed  by  taking  the  cross-product  of 
two  tangent  vectors.  The  tangent  vectors  can  be  obtained  by  differentiation  of  the 
parametric  form  of  the  equation  of  the  surface. 


5.  Tesselation  of  the  Gaussian  sphere:  Orientation  Histograms 

It  is  useful  to  divide  the  sphere  up  into  cells  in  order  to  represent  the  information 
on  the  Gaussian  sphere  in  a  computer.  Ideally  the  cells  should  satisfy  the  following 
criteria: 

1.  All  cells  should  have  the  same  area. 

2.  All  cells  should  have  the  same  shape. 

3.  The  cells  should  have  regular  shapes  that  are  compact. 

4.  The  division  should  be  fine  enough  to  provide  good  angular  resolution. 

5.  For  some  rotations  the  cells  should  be  brought  into  coincidence  with  themselves. 

Cells  which  are  compact  combine  information  only  from  surface  patches  which  have 
nearly  the  same  orientation.  Elongated  cells  of  the  same  area  combine  information 
from  surface  patches  which  have  more  widely  differing  orientations.  The  area  of  a 
regular  polygon  with  n  sides  inscribed  in  a  circle  of  radius  r  is 

2[  sin(27r/n)' 

(27r/n)  J‘ 

So  the  area  of  a  hexagon  incribed  in  a  circle  is  (3\/3/2)r2,  twice  that  of  a  triangle 
inscribed  in  the  same  circle.  Tesselations  with  near-triangular  cells  will  thus  combine 
information  from  orientations  which  are  \fl  times  as  far  from  the  average  as  do 
tesselation  using  near-hexagonal  cells. 
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Figure  15.  The  Gaussian  sphere  can  be  divided  into  cells  along  meridians  and 
es  of  longitude.  The  resulting  cells  do  not  have  the  same  areas  however  and  only 
align  with  each  other  for  certain  rotations  about  the  axis  through  the  poles. 

If  cells  occur  in  a  regular  pattern,  the  relationship  of  a  cell  to  its  neighbors  will 
be  the  same  for  all  cells.  Such  arrangements  are  to  be  preferred.  Unfortunately,  it 
is  not  possible  then  to  simultaneously  satisfy  the  criteria  listed  above. 

A  simple  tesselation  consists  of  a  division  into  latitude  bands,  each  of  which  is 
then  further  divided  along  longitudinal  strips  (Figure  15).  The  cells  could  be  made 
more  nearly  equal  in  area  by  having  fewer  at  higher  latitudes,  or  by  making  the 
latitude  bands  wider  there,  or  both.  One  advantage  of  this  scheme  is  that  it  makes 
it  easy  to  compute  to  which  cell  a  particular  surface  normal  should  be  assigned. 
Still,  this  arrangement  does  not  come  close  to  satisfying  the  criteria  stated  above. 
In  particular,  the  cells  are  brought  into  alignment  only  for  a  few  rotations  about 
the  axis  of  the  globe.  Rotations  about  any  other  axis  can  not  bring  the  cells  into 
alignment. 

5.1.  Tcsselations  Based  on  Regular  Polyhedra 

Better  tesselations  may  be  found  by  projecting  regular  polyhedra  onto  the  unit 
sphere  after  bringing  their  center  to  the  center  of  the  sphere  [31]  Regular  polyhedra 
are  uniform  and  have  faces  which  are  all  of  one  kind  of  regular  polygon  (They  are 
also  called  the  Platonic  solids).  [8,  10,  19,  23,  25,  26,  30].  The  vertices  of  a  regular 
polyhedron  arc  congruent.  A  division  obtained  by  projecting  a  regular  polyhedron 
has  the  desirable  property  that  the  resulting  cells  all  have  the  same  shape  and  area. 
Also,  all  cells  have  the  same  geometric  relationship  to  their  neighbors.  In  the  case  of 


the  dodecahedron,  the  cells  are  even  fairly  well  rounded.  The  dodecahedron,  however, 
has  only  twelve  cells  (Figure  16a).  Even  the  icosahedron,  with  twenty  triangular 
cells,  provides  too  coarse  a  sampling  of  orientations  (Figure  16b).  Furthermore, 
its  cells  are  not  well  rounded.  Unfortunately  there  are  only  five  regular  solids 
(tetrahedron,  hexahedron,  octahedron,  dodecahedron,  and  icosahedron). 

One  can  go  a  little  further  by  considering  semi-regular  polyhedra.  A  semi-regular 
polyhedron  has  regular  polygons  as  faces,  but  the  faces  are  not  all  of  the  same 
kind  [8,  10,  19,  23,  25,  26,  30]  (They  are  also  called  the  Archimedean  polyhedra). 
As  for  regular  polyhedra,  the  vertices  are  congruent.  There  can  be  either  two  or 
three  different  types  of  faces  and  these  have  different  areas.  An  illustration  of  a 
tessclation  using  a  semi-regular  polyhedron  is  provided  by  a  soccer  ball  (Figure 
17a).  it  is  based  on  the  truncated  icosahedron,  a  semi-regular  polyhedron  which 
has  12  pentagonal  faces  and  20  hexagonal  faces. 

Unfortunately  there  are  only  13  semi-regular  polyhedra  (The  five  truncated 
regular  polyhedra,  cuboctahedron,  icosidodecahedron,  snub  cuboctahedron,  snub 
icosidodecahedron,  truncated  cuboctahedron,  rhombicuboctahedron,  truncated 
icosidodecahedron,  and  the  rhombicosidodecahedron).  Overall,  these  objects  do 
not  provide  us  with  fine  enough  tesselations,  The  snub  icosadodecahedron  has  the 
largest  number  of  faces,  but  each  of  its  80  triangles  is  much  smaller  than  each  of 
its  12  pentagons. 

The  edges  of  a  semi-regular  polyhderon  are  all  the  same  length.  One  consequence 
of  this  is  that  the  different  types  of  faces  have  different  areas.  The  area  of  a  regular 
polygon  of  n  sides  and  edge-length  e  equals 
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Figure  17.  Tesselation  of  the  Gaussian  sphere  using  (a)  the  truncated  icosahedron 
and  (b)  the  pentakis  dodecahedron. 


4  tan(7r/n)’ 


so  it  is  very  roughly  proportional  to  n2.  This  is  a  problem  generally  with  semi-regular 
polyhedra.  It  is  sometimes  possible  to  derive  a  new  polyhedron  which  has  the  same 
adjacency  relationships  between  faces  as  a  given  semi-regular  polyhedron  but  also 
has  faces  of  equal  area.  The  shapes  of  some  of  these  faces  then  are  no  longer  regular, 
however. 

if  we  desire  a  finer  subdivision  still,  we  can  consider  splitting  each  face  of 
a  given  tesselation  further  into  triangular  facets.  If,  for  example,  we  split  each 
pentagonal  face  of  a  dodecahedron  into  five  equal  triangles  we  obtain  a  pentakis 
dodecahedron  with  60  faces  (Figure  17b).  This  happens  to  be  the  dual  of  the 
truncated  icosahedron,  discussed  above.  If  we  apply  this  method  instead  to  the 
truncated  icosahedron  we  construct  an  object  with  180  faces.  This  object,  as  well 
as  the  pentakis  dodecahedron,  form  suitable  bases  for  further  subdivision,  as  we 
shall  show  later. 

To  see  how  fine  a  division  we  might  need,  let  us  calculate  the  angular  spread 
of  surface  normals  which  map  into  a  particular  cell.  If  there  are  n  equal  cells,  then 
each  one  will  have  area 

/I  =  (4  7r)/ n, 

since  the  total  area  of  the  unit  sphere  is  47r  (This  area  equals  the  solid  angle  of  the 
cone  formed  by  the  cell  when  connected  to  the  center  of  the  sphere).  The  shape 
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which  minimizes  the  angular  spread  for  given  surface  area  is  the  circular  disc.  The 
area  of  a  circular  disc  on  the  unit  sphere  is 

A  =  27r(  1  —  cos0), 

where  0  is  the  half-angle  of  the  cone  formed  by  the  disc  when  connected  to  the 
center  of  the  sphere.  I1'  0  is  small,  the  area  can  be  approximated  by 

A  %  7 \0l . 

Thus  if  there  are  many  cells  and  if  they  could  be  made  circular,  the  angular  spread 
would  be 

6^  2/yfi. 

The  best  we  can  hope  for,  however,  are  near-hexagonal  cells.  The  area  of  a  hexagon 
inscribed  in  a  circle  of  radius  r  is  (3\/3/2 )r2,  as  already  mentioned.  The  area  of  the 
circle,  rrr*.  is  about  20  °c  more.  So  a  hexagonal  shape  has  a  spread  which  is 

2T  =  1.0996. . . 

\  3\/3 

as  large  as  that  of  a  circular  shape  of  equal  area.  A  lower  bound  on  the  angular 
spread  for  a  tessellation  with  n  cells  then  is 

|^4~_  2.1993... 

\  3^3 n  \fn 

For  n  —  GO,  for  example,  the  spread  is  at  greater  than  16.2  degrees.  One  should 
also  remember  that  the  spread  for  triangular  cells  is  even  more,  namely  >/2  times 
that  for  hexagonal  cells. 

5.2.  Geodesic  Domes 

To  proceed  further,  we  can  divide  the  triangular  cells  into  four  smaller  triangles 
according  to  the  well  known  geodesic  dome  constructions  [19.  27,  31].  We  attain  high 
resolution  by  relenting  on  several  of  the  criteria  given  above  (Figure  18).  Specifically, 
the  colls  of  a  geodesic  tessolation  do  not  all  have  the  same  area  and  shape.  The  cells 
are  also  not  compact,  being  shaped  like  (irregular)  triangles.  The  duals  of  geodesic 
domes  are  better  in  tins  respect,  since  they  have  facets  that  are  mostly  (irregular) 
hexagons,  with  a  dozen  (regular)  pentagons  thrown  in.  Tesselations  of  arbitrary 
fineness  can  be  constructed  in  this  fashion.  The  pentakis  dodecahedron  is  a  good 
starting  point  for  a  geodesic  division,  as  is  the  object  constructed  earlier  from  the 
truncated  icosahedron  by  dividing  the  faces  into  triangles. 

Fvach  of  the  edges  of  the  triangular  cells  of  the  original  polyhedron  are  divided 
into  /  sections,  where  /  is  called  the  frequency  of  the  geodesic  division.  The  result 
is  that  each  face  is  divided  into  f2  (irregular)  triangles.  Tesselations  where  the 
frequency  is  a  power  of  two  are  particularly  well  suited  to  the  method  suggested 
here,  as  we  sec  next. 

One  has  to  be  able  to  efficiently  compute  to  which  cell  a  particular  surface 
normal  belongs.  In  the  case  of  the  tesselations  derived  from  regular  polyhedra,  one 
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Figure  18.  Tesselation  of  the  Gaussian  sphere  using  a  frequency  four  geodesic 
tesselation  based  on  the  icosahedron  (There  are  16  X  20  =  320  faces). 


Prst  computes  the  dot-product  of  the  given  unit  vector  and  the  vector  to  the  center 
of  each  ceil  (These  reference  vectors  correspond  to  the  vertices  of  the  dual  of  the 
original  regular  polyhedron).  This  gives  one  the  cosine  of  the  angle  between  the 
two.  The  closest  reference  vector  is  the  one  which  gives  the  largest  dot-product. 
The  given  vector  is  then  assigned  to  the  cell  corresponding  to  that  reference  vector. 

In  the  case  of  a  geodesic  dome,  it  is  possible  to  proceed  hierarchically, 
particularly  if  the  frequency  is  a  power  of  two.  The  geodesic  dome  is  based  on  some 
regular  polyhedron.  The  appropriate  facet  of  this  polyhedron  is  found  as  above. 
Next,  one  determines  into  which  of  the  triangles  of  the  first  division  of  this  facet 
the  given  unit  normal  falls.  This  can  be  done  by  considering  which  dot-product  has 
the  second  largest  value.  No  new  dot-products  need  to  be  computed.  The  process 
is  then  repeated  with  the  four  triangles  into  which  this  facet  is  divided,  and  so  on. 
In  practice,  lookup  table  methods  can  be  used,  which,  while  not  exact,  are  very 
quick. 

Let  the  area  occupied  by  one  of  the  cells  on  the  Gaussian  sphere  be  w  (in 
the  case  of  the  icosahedron  w  =  47r/20).  The  expected  number  of  surface  normals 
mapped  into  a  cell  equals 

pw\G\, 

for  a  convex  object,  where  is  the  average  of  6'(£,  r;)  over  the  cell. 
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It  is  clear  that  the  extended  Gaussian  image  can  be  computed  locally.  One 
simply  counts  the  number  of  surface  normals  that  belong  in  each  cell.  The 
expression  for  the  Gaussian  curvature,  on  the  other  hand,  includes  first  and  second 
partial  derivatives  of  the  surface  function.  In  practice,  estimates  of  derivatives  are 
unreliable,  because  of  noise.  It  is  important  therefore  that  the  extended  Gaussian 
image  can  be  computed  without  estimating  the  derivatives. 

The  values  in  the  cells  can  be  thought  of  as  an  orientation  histogram.  It  has 
recently  be  brought  to  rny  attention  that  this  is  analogous  to  a  scheme  used  for 
histogramming  directions  of  dendrites  on  neurons  [5]. 

The  result  can  be  displayed  graphically  using  normal  vectors  on  each  of  the 
cells  to  represent  the  weight  of  the  accumulated  surface  normals.  A  frequency 
two  sub-division  of  the  pentakis  dodecahedron  provides  240  cells,  enough  for  most 
practical  purposes  (Figure  19).  The  angular  spread  in  this  case  is  about  1 1 .5  degrees. 
An  alternative  way  to  present  the  extended  Gaussian  image  graphically  is  by  means 
of  a  grey-level  image  where  brightness  in  each  cell  is  proportional  to  the  count. 
The  surface  of  the  Gaussian  sphere  may  be  projected  stereographically  instead  of 
orthographically  in  order  to  preserve  the  shapes  of  the  cells.  Their  areas  will  be 
distorted  however. 

6.  Solids  of  Revolution 

In  the  case  of  the  surface  of  a  solid  of  revolution,  the  Gaussian  curvature  is 
rather  easy  to  determine.  The  solid  of  revolution  can  be  produced  by  rotating  a 
(planar)  generating  curve  about  an  axis  (Figure  20).  Let  the  generating  curve  be 
specified  by  the  perpendicular  distance  from  the  axis,  r(s),  given  as  a  function  of 
arc  length  s  along  the  curve.  Let  0  be  the  angle  of  rotation  around  the  axis.  Now 
consider  the  Gaussian  sphere  positioned  so  that  its  axis  is  aligned  with  the  axis 
of  the  solid  of  revolution.  Let  £  be  the  longitude  and  tj  be  the  latitude  on  the 
Gaussian  sphere. 

We  can  let  £  correspond  to  6.  That  is,  a  point  on  the  object  produced  when 
the  generating  curve  has  rotated  through  an  angle  9  has  a  surface  normal  that  lies 
on  the  Gaussian  sphere  at  a  point  with  longitude  £  =  9. 

6.1.  Gaussian  Curvature  of  Solid  of  Revolution 

Consider  a  small  patch  on  the  Gaussian  sphere  lying  between  £  and  £  -(-  6£  in 
longitude  and  between  rj  and  rj  6 rj  in  latitude.  Its  area  is 

cos  r?  StSrj. 

We  need  only  determine  the  area  of  the  corresponding  patch  on  the  object.  It  is 

r606s, 

where  6s  is  the  change  in  arc  distance  along  the  generating  curve  corresponding 
to  the  change  6r]  in  surface  orientation.  The  Gaussian  curvature  is  the  limit  of  the 
ratio  of  the  two  areas  as  they  tend  to  zero.  That  is, 
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Figure  19.  Orientation  Histogram  collected  on  a  geodesic  dome  derived  from 
the  pentakis  dodecahedron  (There  are  12  X  5  X  4  =  240  faces).  This  is  a  discrete 
approximation  of  the  extended  Gaussian  image.  1  he  length  of  the  vector  attached 
to  the  center  of  a  cell  is  proportional  to  the  number  of  surface  normals  on  the 
surface  of  the  original  object  which  have  orientations  falling  within  the  range  of 
directions  spanned  by  that  cell. 


K  = 


lim 

A  *1  — •  0 


cos  rj  6£  6rj 
t  66  6s 


..  COS  T)  6ti 

lim  - — 

•Sr/— ►  o  r  6s 


cos  rj  dr) 
r  ds  ’ 


since  6£  =  69.  The  curvature  of  the  generating  curve,  kq,  is  just  the  rate  of  change 
of  direction  with  arc  length  along  it  [11,  20,  24].  So 


and  hence 


K  _  *C  COST) 

r 

It  is  easy  to  see  that  (Figure  21),  sin  r)  =  — rs,  where  r,  is  the  partial  derivative  of 
r  with  respect  to  s.  Differentiating  with  respect  to  s  we  get 


26 


Horn 


EG! 


Figure  20.  A  solid  of  revolution  can  be  generated  by  rotating  a  curve  around  an 
axis.  The  curve  can  be  specified  by  giving  the  distance  form  the  axis  as  a  function 
of  the  arc  length  along  the  curve. 


dr\  d  . 
cm"  *  = 

and  so  we  obtain  the  simple  formula 


In  the  case  of  a  sphere  of  radius  R,  for  example,  we  have  r  =  Rcos(s/R)  for 
—{■n/2)R  <  s  <  +(7r/2 )R.  Thus  raa  =  —  (r//?2)  and  K  =  l/R2. 

For  some  purposes  it  is  more  useful  to  express  the  radius  r  as  a  function  of 
the  distance  along  the  axis,  rather  than  as  a  function  of  arc  length  along  the  curve. 
Let  the  distance  along  the  axis  be  denoted  by  z.  It  is  easy  to  see  that  (Figure  21), 
tan  rj  =  —rz,  and  so,  differentiating  with  respect  to  s, 
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Figure  21.  This  figure  shows  the  relationships  between  infinitesimal  increments  in 
arc  length  along  the  curve,  distance  from  the  axis  of  rotation  and  distance  along 
the  axis  of  rotation. 


2  drf  d  dz 

sec  IT  ~  T"( — r*)  =  —r**T' 
as  as  ds 

where  from  the  figure  we  see  that  cos  r/  =  zs,  so  that 


dr/  4 

cos  V  =  —  cos  7]  =  —rzz  cos  rj, 
ds  ' 


and  finally 


K 


r(l  +rir 


since 


sec2  r)  =  1  -j-  r2. 


6.2.  Alternate  Derivation  of  Gaussian  Curvature  of  a  Solid  of  Revolution  (*) 

We  first  need  to  review  Meusnier’s  theorem  [11,  20,  24].  Consider  a  normal 
section  of  a  surface  at  a  particular  point.  It  is  obtained  by  cutting  the  surface 
with  one  of  the  planes  including  the  local  normal  (Figure  22).  Suppose  that  the 
curvature  of  the  curve  in  which  the  surface  cuts  this  plane  is  Now  imagine 
tilting  the  plane  away  from  the  normal  by  an  angle  rj  (using  the  local  tangent  as 
an  axis  to  rotate  about).  The  new  plane  will  cut  the  surface  in  a  curve  with  higher 
curvature.  In  fact,  it  can  be  shown  that  the  new  curve  has  curvature 
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It  is  easy  to  see  this  in  the  case  of  a  sphere,  since  a  plane  including  the  center  cuts 
the  sphere  in  a  great  circle,  while  an  inclined  plane  cuts  it  in  a  small  circle  of  radius 
proportional  to  the  cosine  of  the  angle  of  inclination. 

Now,  let  us  return  to  the  surface  of  revolution.  It  is  not  hard  to  show  that 
one  of  the  principal  curvatures  at  a  point  on  the  surface  will  correspond  to  a  cut 
through  the  surface  by  a  plane  which  includes  the  axis  of  revolution.  The  curve 
obtained  in  this  way  is  just  the  generating  curve  of  the  solid  of  revolution.  So  one 
of  the  two  principal  curvatures  is  equal  to  the  curvature  kq  of  the  generating  curve 
at  the  corresponding  point. 
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Now  consider  a  plane  perpendicular  to  the  axis  of  revolution  through  the 
same  surface  point  (Figure  23).  It  cuts  the  surface  in  a  circle.  The  curvature  in 
this  plane  equals  (1  /r),  where  r  is  the  radius  of  the  solid  of  revolution  at  that 
point.  This  horizontal  plane,  however,  is  not  a  normal  section.  Suppose  that  the 
normal  makes  an  angle  t)  relative  to  this  plane  (The  local  tangent  plane  also  makes 
an  angle  f?  relative  to  the  axis  of  revolution).  Now  construct  the  plane  including 
the  local  normal  which  intersects  the  horizontal  plane  in  a  line  perpendicular  to 
the  axis.  This  plane  will  be  inclined  rj  relative  to  the  one  we  have  just  studied. 
It  also  produces  the  second  principal  normal  section  sought  after.  By  Meusnicr’s 
theorem  we  see  that  the  curvature  of  the  curve  found  in  this  normal  section  is 
Kjs/  =  (l/r)cosr;.  Finally,  the  Gaussian  curvature  is  found  by  multiplication  to  be 

K  _  cos  T) 
r 

In  the  case  of  a  sphere  of  radius  R,  for  example,  we  have  r  =  R  cos  rj  and  kq  —  1  / R, 
so  that  K  —  l//?2,  as  expected. 

To  make  this  result  more  usable,  erect  a  coordinate  system  with  the  z-axis 
aligned  with  the  axis  of  revolution.  The  generating  curve  is  given  as  r(z).  Let 
the  first  and  second  derivatives  of  r  with  respect  to  z  be  denoted  by  r2  and  r„ 
respectively.  It  is  easy  to  see  that  (Figure  21),  tan  r]  =  rz,  so  that 

1 

cos  rj  =  - . 

yjl+'l 


Furthermore, 


kg  — 


(1  +  r2)3/2’ 


so  that  finally, 


K  = 


In  order  to  use  this  result  in  deriving  extended  Gaussian  images  it  is  necessary  to 
identify  points  on  the  surface  with  points  on  the  Gaussian  sphere.  Suppose  that  we 
introduce  a  polar  angle  9  such  that 


x  =  r  cos  9  and  y  =  rsin0. 


Then  a  unit  normal  to  the  surface  is  given  by 


(cos  9 ,  sin  9,  —r,)T 

\A  +r2 

Equating  this  to  the  unit  normal  on  the  Gaussian  sphere, 

(cos£  cos  rj,  sin  $  cos  r/,  sin  rj)T , 


we  get 


£  —  6  and  tan  rj  =  —rt. 
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6.3.  Extended  Gaussian  Image  of  a  Torus 


As  an  illustration  we  will  now  determine  the  extended  Gaussian  image  of  a 
torus.  Let  the  torus  have  major  axis  R  and  minor  axis  p  (Figure  24).  A  point  on 
the  surface  can  be  identified  by  6  and  s,  where  6  is  the  angle  around  the  axis  of 
the  torus,  while  s  is  the  arc  length  along  the  surface  measured  from  the  plane  of 
symmetry.  Then, 

r  =  R  + pcos(a/p), 

and 

r„  =  — (l/p)cos(s/p), 

so  that 
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Figure  24.  A  torus  obtained  by  spinning  a  circle  around  an  axis.  The  resulting 
object  is  not  convex.  Its  extended  Gaussian  image  can  be  computed  nevertheless. 


K  _  ___  1  cos [s/p) 

r  p  R  4-  pcos(s/p)' 

Two  points,  P  and  P'  (Figure  24),  separated  by  n  in  9,  have  the  same  surface 
orientation  on  the  torus.  The  surface  normal  at  one  of  these  places  points  away 
from  the  axis  of  rotation,  while  it  points  towards  the  axis  at  the  other  place. 
Accordingly,  two  points  on  the  object, 

{9,s)  =  {Z,PV)  and  {9,  s)  =  (£  +  tt,  p  (tt  —  rj)), 

correspond  to  the  point  (£,t?)  on  the  Gaussian  sphere.  The  curvatures  at  these  two 
points  have  opposite  signs 


K. p  — 


1  COST? 


p  R  -f-  P  COS  T? 


and 


1  COST? 
p  R  —  p  COS  T? 


The  torus  is  not  a  convex  object,  so  more  than  one  point  on  its  surface  contributes 
to  a  given  point  on  the  extended  Gaussian  image.  If  we  add  the  absolute  values  of 
the  inverses  of  the  curvature  we  get 


1 


1 


If  we  had  added  the  inverses  algebraically  instead  we  would  have  obtained 

iT+  +  k1  =  21,1  • 

which  is  twice  the  result  for  a  sphere  of  radius  p. 

The  same  results  could  have  been  found  using 

K  —  KGcosr 1 
r  ' 

since  kq  —  —\/p  and  r  =  R  ^  pcosi),  so  that 

±COST] _ 

p{R±pcosr))' 

The  extended  Gaussian  image  of  a  torus  has  singularities  at  the  poles.  These 
correspond  to  the  two  rings  on  which  the  torus  would  rest  if  it  were  dropped  onto 
a  plane.  All  of  'he  points  on  one  of  these  rings  have  the  same  surface  orientation. 

We  can  think  of  the  Gaussian  sphere  as  covered  by  two  stieets  of  a  Riernann 
surface,  one  corresponding  to  the  inner  half  of  the  torus,  closer  to  its  axis  of 
symmetry,  the  other  corresponding  to  the  outer  haif.  The  two  sheets  arc  connected 
to  one  another  at  the  poles,  branch  points  corresponding  to  the  two  rings  mentioned 
above.  There  the  Gaussian  curvature  changes  sign. 

We  may  also  note  at  this  point  that  all  tori  with  the  same  surface  area,  (47r2p/?), 
have  the  same  extended  Gaussian  image. 

6.4.  The  (Tuque  Convex  Object  w  ith  G(£,r?)  —  2  sec  77  (*) 

While  all  tori  with  surface  area  4tt2  have  the  same  extended  Gaussian  image, 

G{i,r])  =  2  seer?, 

there  is  only  one  convex  object  which  has  that  extended  Gaussian  image.  It  is  a 
solid  of  revolution  since  G(^,r?)  is  independent  of  £.  So  we  have  on  the  one  hand 

K  =  1/2  cos  r?, 


and  on  the  other  hand 


so  that 


K  _  COST? 


=  r/2. 

The  equation  states  that  the  curvature  of  the  generating  curves  varies  linearly  with 
the  distance  from  the  axis  of  rotation.  This  deceptively  simple  equation  represents 
a  non-linear  second  order  differential  equation  for  r  in  terms  of  z  since 


KG  =  — 


so  that 
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Now 

•  d  1  _  T*r j» 

d*s/T+7i  (i  +  i)3/! 

and 

d  r2  _  rrz 
dz  4  2 

so  that 

d  1  dr2 

dz  n  ,  ^  dz  4 

v1  +  r* 

or 

1 _  =  r2+c2 

\A+r*  4 

where  c2  is  a  constant  of  integration.  We  now  have  reduced  the  problem  to  a 
non-linear  first  order  differential  equation  for  r  in  terms  of  z.  If  the  object  is  to  be 
convex  and  smooth  at  its  poles,  we  expect  r2  — ►  c©  as  r  ->  0.  Thus  c  =  0.  Next 
note  that  the  term  on  the  left  equals  cost/.  So  we  also  have 

r2 

cos  n  =  — 

4 

or,  using  an  earlier  expression  for  k.q, 

*G  =  —y/cosrj. 

This  is  an  implicit  equation  for  the  curve  of  least  energy  [14]!  The  curve  of  least 
energy  is  the  curve  which  minimizes  the  integral  of  the  square  of  the  curvature  kq. 
It  can  be  solved  for  z  in  terms  of  r  to  yield 

z  =  Vr2[2E(cos~1(r/2),  l/\/2)  —  F(cos—1(r/2),  l/\/2)). 

where  E  and  F  are  incomplete  elliptic  integrals.  If  wc  let  s  be  the  arc  length  along 
the  curve  wc  can  also  write  the  solution  in  Whewell  form 

S  =  x/2F{cOS_1  ^/COS  T],  l/y/2), 

or  Cesaro  form 

s  =  x/^Ffcos-^ — tic),  l/y/2). 

The  length  of  the  curve  from  the  pole  to  the  equator  is 

y/2K  [l/y/2)  =  /A))  = 

2v2tt 

where  K  is  the  complete  elliptic  integral  of  the  first  kind,  and  T  is  the  gamma- 
function  [20].  The  height  from  equator  to  the  pole  is 
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while  the  maximum  radius  is 

H  =  2. 

The  minimum  radius  of  curvature  equals  one,  so  that  a  circle  tangent  at  the 
outermost  point  is  also  tangent  at  the  origin  [14].  This  circle,  when  rotated  about 
the  vertical  axis,  produces  a  torus  with  the  same  extended  Gaussian  image  (Figure 
25).  Both  objects  have  total  surface  area  47r2. 

7.  Gaussian  Curvature  in  the  General  Case 


When  the  object  is  not  a  solid  of  revolution  we  need  to  work  a  little  harder 
to  obtain  the  Gaussian  curvature.  Let  x  =  x(u,  u),  y  =  y(u,  v),  and  z  —  z{u ,  v)  be 


.  / 


Horn 


EGI 


parametric  equations  for  points  on  a  given  surface.  Let  r  =  ( x,y,z)T  be  a  vector 
to  a  point  on  the  surface.  Then 

dr  dr 

—  ~  and  tv  — 
ou  a  v 

are  two  tangents  to  the  surface,  as  already  noted  earlier.  The  cross  product  of  these 
two  vectors, 

n  =  ru  X  r„, 

will  be  perpendicular  to  the  local  tangent  plane  (Figure  M).  The  length  of  this 
normal  vector  squared  equals 

n2  =  n  ■  n  =  (ru  ■  rt4)(rv  ■  rv)  —  (ru  ■  rv)2 

since  (a  X  b)  •  (c  X  d)  =  (a  ■  c)(b  ■  d)  —  (a  •  d)(b  •  c).  A  unit  vector  n  =  n / n  can  be 
computed  using  this  result. 


7.1.  Gaussian  Curvature  from  Variation  in  Normals  (*) 


The  Gaussian  curvature  is  the  limit  of  the  ratio  of  the  area  of  a  patch  on  the 
Gaussian  sphere  to  the  area  of  the  corresponding  patch  on  the  surface,  as  the  area 
shrinks  to  zero.  Consider  an  infinitesimal  triangle  formed  by  the  three  points  on 
the  surface  corresponding  to  (u,r),  (ix Cu,  <.  ),  and  lu,e  -]-£r).  The  lengths  o:' two 
sides  of  this  triangle  are 

jrtii<5u  and  |r„|6v 

while  the  sine  of  the  angle  between  these  sides  equals 


jr„  X  rvj 
M  |r„| 

so  that  an  outward  normal  with  size  equal  to  the  area  of  the  triangle  is  given  by 


-(r„  X  Tv)6u6v  = 


n  6u6v. 

2 


To  v.,  ermine  the  area  of  the  corresponding  triangular  patch  on  the  Gaussian  spheae 
we  need  to  find  the  unit  surface  normals  at  the  three  points.  The  unit  surface 
normals  will  oe 


h,  h  -f  hu  Su,  and  h  4*  Sv 


if  we  ignore  terms  of  higher  order  in  £u  and  fiv.  Here  hu  and  nv  are  the  partial 
derivatives  of  n  with  respect  to  u  and  v.  V  .  that  utJ  and  nt.  are  perpendicular  to 
h.  The  area  of  the  patch  on  the  Gaussian  ">tu-re  equals  the  magnitude  of 


^(nu  X 

by  reasoning  similar  to  that  used  in  determining  the  area  of  the  original  patch  on 
the  given  surface.  We  need  to  find  nu  and  nv  to  compute  this  area.  Now 


d  n  nnu  —  nn„ 
du  n  n 2 
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From  n2  —  n 

■  n  we  get 

n  nu  =  a  •  nu, 

so  that 

(n  •  n)uy  —  (n  •  nu)n 

(n  X  nu)  X  n 

Hu  — 

n3 

n3 

and 

(n  •  n)nt  —  (n  ■  n.  )n 

(n  X  nv)  X  n 

Ht  — 

n3 

n3 

since  (a  X  b) 

X  c  —  (a  • 

c)b  —  (b  •  c)a.  'fhen 

VI 

n„  X  nv  =  ,  |(n  •  n)(nu  X  nv)  -f  (n  ■  nj(n„  X  n)  -f  (n  •  njfn  X  n„)l 

n° 

or 

1  r  , 

n«  X  Dj  -  --  (n  nu  nvjn 
nq 

.since  fab  clp  =  (a  p)(b  X  c)-v(b-p)(c  X  a)  +  (c  p)(aXb),  where  [abrj  =  (aXb]  c. 

This  shows  that  the  patch  on  the  Gaussian  sphere  has  the  same  orientation  as 
the  patch  on  the  surface,  as  it  should.  An  outward  pointing  normal  of  size  equal  to 
the  area  is  given  by 


r  -7  [n  n„  n„]n  SuSv. 

2  n  * 

The  ratio  of  the  two  areas,  the  Gaussian  curvature,  is  then  just 


Now 
so  that 


K  = 


[n  nu  nv] 


n  =  ru  X  rv, 


and 


nu  =  ruu  X  rv  -f  ru  X  r„ 


Hy  -  tyy  X  Fy  —j~  X  Tyv 

Now  using  (a X  b)  x  (c  X d)  =  [a  b  d]c—  [a b  c]d  or  (a X  b)  X  (c  X  d)  =  [ac  d]b  —  [be  dja 
we  get 

nu  X  ny  —  [ruu  rv  ryy]rv  -p  [ruu  rv  ryy]ru  [rUy  rt,  Tyjryy -f-  [ru  TyV  rv]rUy -j-  [ru  ryy  Tyy]ru 
so  that 


[n  Dy  ny]  =  H'(n„X  nv)  =  [r„  rv  ryU][ry  rv  rVy]  —  [r„  rv  ryV]2, 


and  finallv 


K  = 


K  Ty  Tyy][ry  Ty  Tyy]  [ryl"v^1iv] 

jru  X  ryj4 


This  result  can  be  used  to  derive  the  expression  for  curvature  of  a  solid  of  revolution 
in  a  more  rigorous  fashion. 
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7.2.  Fundamental  Forms  of  a  Surface  (*) 
Let,  as  before, 


-  —  L*  x  rv 
|r„  X  r„| ’ 

be  the  unit  surface  normal  vector.  The  first,  fundamental  form  of  a  surface  gives 
the  square  of  the  element  of  distance  as  [20] 

ds 2  =  |dr|2  —  E(u,  v)  du 2  -j-  2F(u,  v)  du  dt  +  G[u,  v)  dv2 . 

The  second  fundamental  form  of  a  surface  gives  the  normal  curvature  using  the 
equation  [20] 


-dr  •  rfii  =  /^u,t')tiu2  -4-  2Af(u,  v)  du  dv  +  ,Y(u,  t>)  dv2. 
'Hie  coelfirents  can  be  expressed  in  terms  of  derivatives  of  r  as  follows, 


E  —  ru 

•ru,  F  =  rurv,  and  G  —  rt, 

•  r„ 

and 

L  —  r^  n,., 

A /  =  ru  •  nv  =  rv  nu,  and  X 

11 

n 

3» 

or, 

1  [t*u  r*  ruu] 

"v/EG-F2’ 

..  ! rti  rt  rut]  ,  •• 

M  m  -  ,  and  A  - 

v/+;  -  f2 

[ru  r v  r„y] 

Veg  -f2 

so  that  [20] 


__  LX  -  M2 

EG  —  F2  ' 

Finally,  if  the  surface  is  given  as  z(x,y),  the  above  reduces  to  the  familiar 

,2 


K  - 

(1  +  2 


'■'ix-yy 

Y 

x 


7.3.  Application  of  the  General  Formula  to  the  Ellipsoid  (*) 

In  the  case  of  the  ellipsoid  we  have,  as  discussed  before, 

r  -  (a  cos  9  cos  o ,  b  sin  9  cos  d,  c  sin  d)7" , 
rg  -----  ( —  a  sin  0  cos  0,  bcos  0  cos  d,  0)7', 
r^  --  (  a  cos  0  sin  o,  —b  sin  0  sin  d,  c  cos  d)7", 

and 

rgg  =  (—a  cos  0  cos  d,  —b  sin  9  cos  d,  0)7", 
rg 0  (a  sin  0sin  d,  —b  cos  0  sin  <j>,  0)7', 

r<pt>  —  (  —  a  cos  0  cos  0,  — b  sin  6  cos  d,  —  c  sin  d)7'- 
A  surface  normal  can  be  found  by  taking  cross  products, 

n  =  rg  X  r$  =  (be  cos  0  cos  d.  ca  sin  9  cos  d,  ab  sin  <j>)T  cos  d, 
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and  the  coefficients  of  t fie  first  fundamental  form  are, 


E  =  r$  •  r<?  =  (a2  sin2  0  +  b~  cos2  0 )  cos2  <p, 

E  —  r$  •  r^>  =  (a2  +  62)  sin  9  cos  0  sin  <6  cos  0, 

G  =  r$  r^,  =  (a2  cos2  0  +  b2  sin2  9)  sin2  4>  +  c‘  cos2  <f>. 


Hence 

EG  —  F2  —  [(6c  cos  9  cos  <j}2  -f  ( ca  sin  9  cos  o)2  -f-  (a6  sin  d>)2]  cos2  <p. 
To  compute  the  coefficients  of  the  second  fundamental  form  we  need, 


[re  r^,  too]  =  n  •  r oo  =  —abc  cos3  <p, 

[re  r$  r$^,j  n  ■  0, 

[re  r#  Tfit)  —  n  ■  r^  =  —abc  cos  <j>, 


and  so 


[re  r^  r^„l[r9  t0  r^J  —  fr#  r^]2  — ■  (abc  cos2  <t>) 


Finally  then, 


and 

LN  -M2 
“  EG'- F T 


LN  -  M2 


(a6c  cos2  <j>)2 
EG  -  F2 


abc 

(be  cos  0cos  <jl>)2  -f-  (ca  sin  0cos  <t>)2  4*  ( a6sin  <p )2 


This  result  was  used  earlier  in  the  discussion  of  the  extended  Gaussian  image  of 
the  ellipsoid. 
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